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ABSTRACT 



We present Very Long Baseline Array (VLBA) H i absorption observations of the core region of the powerful radio galaxy Cygnus A. 
These data show both broad (FWHM = 231 ± 21km s" 1 ) and narrow (FWHM < 30km s" 1 ) velocity width absorption components. 
The broad velocity absorption shows high opacity on the counter-jet, low opacity against the core and no absorption on the jet side. 
We argue that these results are most naturally explained by a circumnuclear H i absorbing disk orientated roughly perpendicular to the 
jet axis. We estimate that the Hi absorbing gas lies at a radius of ~ 80 pc has a scale height of about 20 pc, density n > 10 4 cm" 3 and 
total column density in the range 10 23 - 10 24 cm" 2 . Models in which the Hi absorption is primarily from an atomic or a molecular 
gas phase can both fit our data. Modelling taking into account the effective beam shows that the broad H i absorbing gas component 
does not cover the radio core in Cygnus A and therefore does not contribute to the gas column that blocks our view of the hidden 
quasar nucleus. If however Cygnus A were observed from a different direction, disk gas on ~ lOOpc radius scales would contribute 
significantly to the nuclear column density, implying that in some radio galaxies gas on these scales may contribute to the obscuration 
of the central engine. We argue that the circumnuclear torus in Cygnus A contains too little mass to power the AGN over > 10 7 yr 
but that material in the outer Hi absorbing gas disk can provide a reservoir to fuel the AGN and replenish torus clouds. The second 
narrow Hi absorption component is significantly redshifted (by 186km s" 1 ) with respect to the systemic velocity and probably traces 
infalling gas which will ultimately fuel the source. This component could arise either within a tidal tail structure associated with a 
recent (minor) merger or be associated with an observed infalling giant molecular cloud. 

Key words, galaxies: active - galaxies: elliptical - galaxies: individual (Cygnus A) - galaxies: kinematics and dynamics - galaxies: 
structure - galaxies: ISM 



1. Introduction 

Circumnuclear obscuring tori/disks are an essential component 
of unified schemes of active galactic nuclei (e.g. lAntonuccH 
1 1 993t. iTadhunterl 12008 ) . Recently there has been much progress 
in obtaining direct evidence for such structures especially 
in Seyfert luminosity objects. This evidence includes mod- 
elling of the IR SEDs from AGN heated dust in clumpy tori 
dNenkova et all l2008bl) and dir ect IR interferometric imaging of 
this d ust on 1 pc - 10 pc scales (iJaffe et alll2004tlTristram et all 
2009). 



On larger scale s adaptive optics IR observations of molecu- 
lar hydrogen lines (Hic ks et all [2009) have revealed geometti- 
cally thick gas at radii 30 pc in Seyferts. Additionally millimette 
interferometry also detects molecular g as in emission on scales 
r — 70 pc (e.g. ISchinnerer et"aQ 12000) albeit in more flattened 
disk-like structures. Such outer disk structures may be contin- 
uous with inner obscuring tori and provide both the fuel and a 
conduit for feeding the central engine. The relationship between 
these circumnuclear disks and obscuring tori is however far from 
clear. 

Radio observations provide another means to study the cir- 
cumnuclear gas environment. This can for instance be achieved 
by V LBI observations of maser emission from molecular gas 
CqI 120051) . free-free absorption from ionised gas and ab- 
sorption from atomic gas (Hi). Examples of the use of the 
latter two tracers include observations of NGC 1275 (e.g. 



I Vermeulen et all 1 19941) . Centaurus A (e.g. Sones et all 119961: 
Moreanti et all 120081) Hvdra A (e.g.lTaylo^ll996l) . NGC 4261 
(Ivan Langevelde et alll2000l:lJones et alll200ll) and 1946 + 708 
dPeck & Tavlorl 1200 lb . Because of their high spatial resolu- 
tion such radio observations are especially suitable for studying 
circumnuclear obscuring matter i n powerful narrow line radio 
galaxies which are expected (see ITadhunterl 120081) to be uni- 
fied via orientation with radio-lou d quasars. There is str ong ev- 
idence from X-ray observations (Hardcastle et all [20091) for the 
expected obscuration by large column densities in the former ob- 
jects, however they are usually too distant and faint for optical 
and IR observations to directly observe the circumnuclear gas. 

A prime target for studies of circumnuclear gas in a lumi- 
nous 'hidden quasar' rad io galaxy is the closest Fanaroff-Riley 
dFanaroff & RilevLll974l) type II (FR-II) radio-galaxy Cygnus A. 
Spectropolarmetric observations of this source which revealed 
a hidden BLR in scattered light (Qgl eet all 119971) were a 
major milestone in the general acceptance of the orientation 
unification scheme for powerful radio galaxies and radio-loud 
quasars. Further evidence for shadowing from a central torus 
come s from the bi - cones observed in both o ptical emission 
lines dJackson et all 1 19981) and IR continuum (ITadhunter et all 
1 1999b . The sharpness of the edges of these bi-cones suggest that 
the inner face of any torus occu r s at radii < 50 pc fro m the 
central engine (Tadhunter, 2008). Tadhunt er et all 0003) have 
from optical/IR emission line observations measured a rota- 
tion curve from gas rotating around the bi-cone/radio-jet axis 
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at r » 300 - 1000 pc allowing a central black hole mass to be 
estimated. These observations may trace the outer parts of a cir- 
cumnuclear gas structure which connects with the inner obscur- 
ing torus. The ultimate origin of this mate rial could be related t o 
the merger activity detected in Cygnus A ((Canalizo et alll2003b . 

The central radio core and inner jets of Cygnus A are 
relatively bright from millimetre to centimetre wavelengths 
allowing searches in absorption to constrain circumnuclear 
gas properties on small (< 100 pc) scales. Molecular ab- 
sorption observations so far give ambiguous results with 
only upp er limits or marginal detec ti ons being r e ported 
(see e.g. iBarvainis & Antonuccil Il994t iFuente et all [2000; 



ISalome & Com bes, 2003; Impellizzer i et all 120061) . VLBI ob- 
servations bv iRrichbau met alJ d 19981) have however found ev- 
idence for ionised circumnuclear gas on scales < 20 pc via 
the detectio n of free-free abso r ption toward the counter-jet. 
Additionally IConwav & Blancol (1995) detected broad Hi ab- 
sorption toward the core in VLA observation which were in- 
terpreted in terms of a circumnuclear disk/torus model, with 
the H i absorption either tracing the small atomic fraction of a 
mainly molecular medium or a purely atomic structure. To bet- 
ter constrain the scale and geometry of this Hi absorbing gas 
we have performed high resolution NRAO Very Long Baseline 
Array (VLBA[3) H i absorption observations . A shor t repor t on a 
initial reduction of this data was given by|C onwav (1999), this 
present paper presents a fuller re-analysis of the data. The or- 
ganisation of this paper is as follows, in Sect. 2 we describe the 
observations while Sect. 3 presents the observational results in- 
cluding modelling of opacity profiles along the source. A discus- 
sion of the results is given in Sect. 4 and a summary in Sect. 5. 
At the redshift of Cygnus A (z = 0.056) for cosmologies with 
H = 73 km s _1 1 mas corresponds to approximately 1 pc, a 
scaling which we adopt throughout this paper. All total reces- 
sion velocities quoted are heliocentric (optical definition). 



2. Observations 

Observations were performed on August 31st 1995 using the ten 
stations of the VLBA plus the phased VLA. Two IFs (left and 
right circular polarisation) with a bandwidth of 12.5 MHz and 
256 channels were centred at the freq uency of the previously ob - 
served Hi absorption (~ 1340 MHz. lConwav & Blancol [1 995). 
A number of bright compact sources were observed as fringe 
finders and for bandpass calibration. The data were correlated in 
Socorro, USA. A standard data reduction using AIPS was per- 
formed including fringe fitting, calibration and flagging of the 
data. Cygnus A lies at a relatively low galactic latitude (b = 5°7) 
and the effects of interstellar scattering are significant at VLBI 
resolution. For this reason four antennas which participated only 
in long baselines, namely Brewster, Hancock, Mauna Kea and 
St. Croix showed no fringes to Cygnus A and so were deleted 
from the subsequent data analysis. Initial amplitude calibration 
was accomplished for VLBA antennas using the recorded sys- 
tem temperature values (which took into account the dominant 
contribution to the noise from the lobes of Cygnus A). For the 
phased VLA the calibrator T ant /T sys values were used with a cor- 
rection applied to take into account the system noise contribu- 
tion from Cygnus A. Offsets in the VLA amplitude calibration 
scale were then corrected by comparing the correlated Cygnus A 
flux densities on long baselines to respectively the VLA and Pie 



Town. Initial continuum images were made via iterative phase 
self-calibration/deconvolution starting from a point source. A 
couple of cycles of amplitude and phase self-calibration were 
performed at the very end to obtain a noise limited continuum 
map. 

Extensive experiments were carried out to determine the op- 
timum uv weighting which gave the best combination of sensi- 
tivity and resolution for both continuum and spectral line. Final 
images were made using uniform weighting with robustness fac- 
tor 0.5; giving rise to an almost circular dirty beam with main 
lobe width FWHM m 25 mas. Since the uv point weights also 
took into account the sensitivity of each baseline (by weighting 
by 1/noise variance) the final continuum and spectral line im- 
ages are dominated by the baselines to the VLA. CLEAN im- 
ages were restored with a circular FWHM 25 mas Gaussian; 
however the effective image resolution is less than this because 
of the effects of foreground interstellar scattering. The observed 
effective resolution (i.e. the clean beam convolved with the inter- 
stellar scattering) is 32.7 mas, as determined by measuring the 
FWHM of the continuum profile perpendicular to the jet at the 
core position - implying an interstellar scattering contribution of 
21 mas. 

The resulting continuum image (see Fig. Q} has rms noise 
of 0.43 mJy beam -1 and shows besides the unresolved core 
a jet and counter-jet structure (see Sect. 3). In principal it is 
possible that the counter-jet structure could be an artefact of 
phase self-calibration starting from a point source. To check 
this possibility we re-mapped the data only allowing flux 
density on the jet side to be included in the model in the 
initial cycles, but in all cases emission on the counter-jet side 
remained. The existence of a counter- jet i s also confirmed by 
other observations (e.g. iKrichbaum et all Il998h . In making 
our final continuum and spectral line images we used data 
self-calibrated against continuum models including both jet 
and counter-jet emission. Before making our spectral line cube 
we removed the continuum contribution using the AIPS task 
UVLIN and then imaged using the same weighting as used for 
our continuum image and the same CLEAN restoring beam. In 
order to increase sensitivity when creating our spectral line cube 
the uv data were averaged in frequency to give a final channel 
separation equivalent to 14.0 km s _1 in velocity or a velocity 
resolution of 28.0 km s _1 after Hanning smoothing. The noise 
achieved in the final line cube was cr rms = 2.36 mJy beam -1 per 
channel. 



3. Results 

3.1. Continuum image 

The 1340 MHz continuum image (Fig. Q] t op panel) is very 
simila r to the 1660 MHz image shown by IKrichbaum et al.l 
dl998l) . We clearly detect the unresolved core, the jet and the 
weaker counter-jet (SE of the core). The PA of the jet is 
105°±2°, in agreement wit h VLB I observations at higher fre- 
quen cies (IKrichbaum et all 1 1998b and the kpc-size jet struc- 
ture (ICariili et all 11991 ). The continuum peak has brightness 
0.16 Jy beam 1 and the total continuum flux recovered in our 
observations is 0.50 Jy. 



1 The VLBA is operated by the National Radio Astronomy 
Observatory which is a facility of the National Science Foundation op- 
erated under co-operative agreement by Associated Universities, Inc. 



3.2. H i absorption 

An integrated absorption spectrum over our spectral line cube 
is shown in Fig Q] bottom left. Broad (Av = 456 km s _1 ) 
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Fig. 1: Top panel: Continuum image at 1340 MHz. The lowest contour is at 2 mJy beam -1 with subsequent contours increasing by 
factors of 2. The effective beam FWHM (see Sect. [2jl is indicated in the lower right corner. Bottom left panel: Integrated absorption 
spectrum from the blanked cube with the velocities shifted to the rest frame of Cygnus A (see Sect. 4.1). The dashed line shows the 
two component Gaussian fit. Bottom right panel: Contours show continuum. Grayscale shows the mean opacity over the rest frame 
velocity range -80 to +170 km s _1 . The thick dark line shows the un-blanked region over which the integrated absorption spectra 
(shown in the bottom left panel) is calculated. 



Hi absorption is detected in the velocity range from 16679 
to 17135 km s^'with the peak being located at v = 
17002 km s _1 (z = 0.05667). The integrated spectrum is well 
fitted by two Gaussian components, yielding centroid velocities 
of 16916 + 10 km s 1 (z = 0.05639 + 0.00003) and 16986 ± 
5 km s- 1 (z = 0.05662 ± 0.00002), with FWHM= 231 ± 21 and 
29 ± 10 km s _1 respectively (after correcting the line widths to the 
rest frame of the host galaxy). The second, narrow component 
has a FWHM velocity width similar to our velocity resolution 
and hence we consider this velocity width as an upper limit of 
the true line width. The presence of two components in the inte- 
grated spectrum suggests that the H i absorption consists of two 
different overlapping structures. Despite small differences in the 
flux scale (for continuum and absorption spectrum) the VLBA 
data agree in absorption width, profile shape and radial velocity 
with th e VLA A- and B -array observations o flConwav & Blanca 
(1995). These authors found slightly different Gaussian fits for 
the absorption spectrum, their data however potentially suffers 
from bandpass calibration and continuum subtraction problems 
(since two IFs — only slightly overlapping in velocity — had to 
be used to achieve sufficient spectral resolution and velocity cov- 
erage). Given these issues we conclude that both data sets agree 
within the noise of both observations so that we have recovered 
the full absorption seen by the VLA. The difference seen in flux 
scale between our VLBI and the published VLA observations are 
likely due to inaccuracies in amplitude calibrations of the VLA. 
This instrument, unlike VLBI, did not record the antenna system 



temperatures which are greatly enhanced due to the presence of 
the bright radio lobes of Cygnus A in the primary beam of each 
antenna, complicating the amplitude calibration. 

Inspection of the data cube shows that the absorption is spa- 
tially resolved and is detectable over 95 mas in angle along the 
radio axis (i.e. ~ 3 effective beam FWHMs). Below the second 
contour of the continuum image it is not possible to constrain the 
H i absorption because the background is too weak. The deep- 
est absorption measured in mJy is toward the counter-jet and 
unresolved core but the highest opacities occur for the broad 
velocity component on the counter-jet side (see Fig.[T] bottom 
right). We find no indication of changes in spectral line absorp- 
tion profiles in directions perpendicula r to the jet axis - which 
is as expected given the small jet width (Krichba um et al.L ll998) 
compared to our effective resolution; this means we need to only 
consider the spectral profile as a function of position along the 
jet axis as shown in Fig. [2] In this figure the top panel shows 
the rotated continuum image while the middle panel shows the 
absorbed flux density (contours) and absorbed flux/continuum 
ratio (colours; note: saturated over 0.2) versus velocity and posi- 
tion along the jet. 

The highest contours in the middle panel of Fig.|2]belong to 
the narrow velocity absorption component at recession velocity 
16986 km s _1 seen against the core. At the same velocity a slight 
extension on the jet-side is detectable, consistent with having 
the same r * 0.1 opacity as seen on the core. Because of the 
rapid fall-off of continuum intensity further along the jet and 
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counter-jet further information about the spatial distribution of 
the narrow velocity component is limited. That is to say that the 
underlying continuum is not strong enough to detect the narrow 
velocity absorption component even if the opacity remains r m 
0.1. 

The broad absorption component in Fig. [2j middle panel 
(seen at velocities around 16900 km s~') is detected from the 
core position out to 65 mas along the counter-jet. Over this spa- 
tial range the broad component has absorbed flux densities (con- 
tour levels) which stay almost constant; given the rapid fall off 
in background continuum this corresponds to a rapid increase 
in line-to-continuum ratio (colours). Over the range of position 
from -50 mas to -65 mas spectra taken show the absorption to 
be fiat bottomed and the absorbed flux density comparable with 
the continuum, both implying H i opacities r » 1 . The absorp- 
tion shows an apparent sharp decrease beyond -65 mas, though 
in a region with very weak continuum (we discuss further the 
reality of this decrease in Sect 3.3). The broad absorption also 
apparently has a small but significant opacity (r w 0.03) at the 
position of the core, however, this may be due only to the lim- 
ited spatial resolution causing "leakage" of absorption onto the 
core position. Quantitative estimates of broad line opacity along 
the source are made in Sect. 13. 31 where the above two points are 
addressed. 

Reinforcing the above description of the two velocity com- 
ponents the bottom panel of Fig. [2] shows the line to continuum 
ratio averaged over different velocity ranges versus position. The 
solid line shows average line-to-continuum for velocity ranges 
where the broad component dominates, while the dashed line is 
for a velocity range where the narrow absorption is most domi- 
nant. The broad component shows rapidly increasing absorption 
along the counter-jet reaching a peak velocity averaged line- 
to-continuum ratio of almost 0.6 at -55 mas (with peak ratios 
within the velocity profile at this position in fact reaching up to 
1 and beyond). The second profile (dashed line) for the veloc- 
ity range where the narrow velocity component normally domi- 
nates has approximately the same mean opacity on the core and 
jet. On the counter-jet side the average opacity over this velocity 
range increases, this is however consistent with the narrow ab- 
sorption component having the same opacity as on the core and 
jet-side but with the average line-to-continuum ratio becoming 
dominated by contamination from the high velocity wings of the 
broad velocity component. 

3.3. Modelling spatial variations in broad line opacity 

The increase in line-to-continuum ratio of the broad absorption 
seen on the counter-jet side (see Sect. 3.2 and Fig. |2j implies a 
rapid increase in opacity. As noted in Sect. 3.2, at the position 
of maximum absorption the spectral profiles are flat-bottomed 
and saturated implying large opacity. To more accurately con- 
vert the line-to-continuum ratios versus velocity to peak opacity, 
and therefore estimate Hi column density variations, we have 
carried out a detailed modelling of the data. In this modelling we 
describe both continuum and peak opacity estimates with a set of 
seven point sources separated by 15 mas along the jet/counter- 
jet axis (i.e. just under half an effective beam FWHM width). 
Additionally by taking into account the effective beam width 
this modelling provides a modest super-resolution of the data, 
important because interstellar scattering limits our spatial reso- 
lution (Sect. 2). Specifically, we are interested in the question of 
whether the weak broad absorption apparently seen at the core 
position (see Fig. Q~]bottom right panel) is real or whether it can 
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Fig. 2: Top panel: Radio continuum image of the core region ro- 
tated such that the jet is pointing to the right and the counter- 
jet to the left. Contour levels are the same as in Fig. Q] Middle 
panel: Position-velocity diagram along the radio axis. Grey 
scale shows the line-to-continuum ratio. Over-plotted in con- 
tours is the absorbed flux density in mJy beam -1 . Contour lev- 
els are -15, -10, -7.5, -5.0 (white) and 5.0 mJy beam -1 (black). 
Bottom panel: Velocity averaged line-to-continuum ratio along 
the source for two different velocity ranges. The solid line is 
for a velocity range dominated by the broad velocity component 
(16720 - 16970 km s -1 ), the dashed line for a velocity range 
(16971 - 17012 km s -1 ) centred on the narrow absorption sys- 
tem velocity. 



be explained by the combination of very strong absorption on 
the counter-jet side combined with limited spatial resolution. 

In our modelling we first made estimates for the underlying 
continuum profile for each model point source (referred to as a 
"pixel") where there was detectable absorption. Continuum in- 
tensities were varied at each of these seven pixels until after con- 
volution with the effective restoring beam the model continuum 
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Table 1 : Results of modelling of continuum and broad H i ab- 
sorption profiles along the jet axis (see Sect- 3.3). 



r 


Fc 


F L 




TO 


(mas) 


(mJy) 


(mJy) 






-60 


4.73 


0.23 


0.05 ±0.12 


0.080 ±0.210 


-45 


3.62 


4.06 


1.12 ± 0.15 


> 3.110 


-30 


0.02 


0.00 


0.00 




-15 


22.78 


4.24 


0.19 ±0.02 


0.320 ± 0.040 





139.16 


1.57 


0.01 ± 0.004 


0.016 ± 0.006 


15 


0.00 


0.00 






30 


20.98 


0.00 


0.00 ± 0.03 


± 0.0478 



Column 1 give the distance from the core of the fitted pixel, neg- 
ative offsets correspond to positions along the counter-jet, posi- 
tive to positions on the jet side. Columns 2 to 4 give respectively 
the fitted continuum flux density, H i absorbed flux density aver- 
aged over fitted frequency range and the line to continuum ratio. 
Column 5 gives the modelled opacity at line centre. 



profile versus position fitted the observations. In a similar way 
we estimated for each pixel the absorbed line area (mJy km s _1 ) 
at velocities in the range 16720 - 16970 km s _1 (where the broad 
absorption component dominates). Pixel values were again ad- 
justed such that after convolution by the effective beam the 
model fitted the observations. Finally an estimate was made 
at each pixel of the peak Hi opacity versus velocity by com- 
bining the pixel-based continuum and spectral line absorption 
estimates. In estimating this quantity we assumed the broad- 
component opacity spectrum was Gaussian with fixed velocity 
centroid and fixed FWHM. 

The resulting fits reproduce (see Fig. |3), to first order, both 
the observed continuum and the absorbed flux profiles along 
the source. Deviations of the fits from the data are likely due 
to the background co ntinuum being more complicated (see 
iKrich baum et al., 1998) than our parameterisation of it. Table. 1 
gives our results and the line centre peak opacity as plotted in 
Fig. [4] Note that no opacity estimate is plotted at positions -30 
and +15 mas because the fitted continuum intensity at these po- 
sitions was zero. We estimate the error of the line-to-continuum 
ratio (Tab. 1, column 4) assuming it is dominated by the rms 
noise in the line data. At position r - -45 mas the lower limit 
on line centre opacity is set by subtracting 2cr from the line-to- 
continuum ratio. 

The peak opacity (column 5) takes into account that part of 
the absorption spectrum is flat bottomed and the attached error 
is calculated based on a lcr rms uncertainty in our line data. Our 
results are consistent with no broad H i absorption on the jet- 
side and any opacity against the core being very low (0.016). 
Moving outward from the core along the counter-jet, the opacity 
increases and peaks at 45 mas (A^ : = 1.6 • 10 21 T sp ; n ) before it 
sharply decreases to 0.078 at 60 mas from the core. Based on 
uncertainty estimates of the absorbed flux this drop in opacity is 
real and not an observational artefact due to insufficient signal- 
to-noise. 



4. Discussion 

4.1. A disk geometry for the broad absorption component 

The fact that the broad velocity width absorption is seen at high 
opacity only against the counter-jet and not against the jet is sim- 
ilar to the situation found for other powerful cores of radio galax- 
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Fig. 3: Modelling of continuum and broad absorption along the 
jet axis (see Sect. 3.3). Top panel: Comparison of model fit and 
data for the observed continuum profile found by adjusting pixel 
values while assuming an effective Gaussian beam of FWHM 
of 32.7 mas. Bottom panel: Comparison of model and data for 
average absorbed line flux over the broad component velocity 
range from 16720 - 16970 km s _1 . In both plots negative offsets 
correspond to counter-jet, positive offsets to jet side. 
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Fig. 4: Estimates of broad absorption opacity at line centre ver- 
sus position along the radio axis after taking into account beam 
effects (see Sect. 3.3). The opacities and errors plotted are taken 
from Tab. 1 . 
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ies observe d in H i absorption, for insta nce in the FR-I source 
NGC 4261 (Ivan Langevelde et all feOQOh . The most natural ex- 
planation for these observations is that the H i absorption lies in 
a flattened structure roughly in a plane perpendicular to the radio 
axis (i.e. a circumnuclear disk). Alternative explanations involv- 
ing foreground clouds in the ISM of the host galaxy seem un- 
likely, requiring within elliptical hosts a large covering factor of 
lOpc-lOOpc sized clouds with uncharacteristically large internal 
velocity dispersion (of order 300 km s _1 ). 

Converting our observations of the location and distribution 
of H i absorption along the counter-jet to an exact of radius in 
the disk depends on the precise orientation of the disk. Over the 
range of plausible disk orientations however this radius varies 
over a fairly narrow range as can be seen by examining two lim- 
its, first when the disk is close to edge-on and when at its max- 
imum plausible tilt. In the former case the H i profile width in 
Fig. [4] measures the scale height of the disk. Sin ce we know the 
central black hole mass (Mbh = 2.5 • 10 9 M^. lTadhunteret all 
I2003h the radius at which the H i absorption occurs can be cal- 
culated (following lKrolik & Begelmanlll988l) via: 

Ah Av 

7-r (1) 

r v rot (r) 

where Ah is the disk scale height, r the radius from the core, 
Av = 98.5 km s _I the velocity dispersion (estimated from the 
Gaussian fit in Sect. 3.3) and v lot (r) the Keplerian rotation veloc- 
ity at radius r from the central black hole. Assuming Ah = 25 pc 
and solving for r we obtain r — 88 pc. In the other limit we 
consider a maximally inclined disk. Based on VLBI observa- 
tio ns of the jet to counter-je t brightness ratio of the parsec scale 
jet dKrichbaum et all 1 1 9981) the jet axis is orientated at an angle 
9 > 80° to the l ine of sight. Based on measured misalignments 
of larger scales (Tadhun ter et alll2003l) we estimate a maximum 
misalignment between disk and radio axes of 30°- which gives 
a maximally inclined disk which is 40° from edge-on. The peak 
H i opacity is observed to occur at projected distance 45 pc along 
the counter-jet which when de-projected corresponds to a radius 
of 70 pc. Using again eq. 1 this gives a scale height of ~ 18 pc. 
Based on the above two limits we will assume in the following 
that the H i peaks at a radius 80 pc from the black hole and has 
a scale height of about 20 pc. The resulting opening angle of the 
circumnuclear H i disk is ~ 14° wh ich is similar to what is foun d 
in other sources (e.g. in NGC 4261 van Lan gevelde et all 12000). 

An important question in considering the feasibility of the 
disk hypothesis is the velocity of the broad absorption relative to 
the systemic velocity of the galactic nucleus. For a disk which 
is perfectly normal to the jet axis these velocities will be the 
same. Our best estimate of the systemic velocity we take to be 
Z = 0.05600 ± 0.00008 (= 16800 ± 24 km s" 1 ) which is the mean 
of six published optical/I R emission line estimates (see Table. 1 
in lTadhunter et all l2003). The centroid of the broad component 
at 16916 km s _1 is therefore 116 km s _1 beyond the mean sys- 
temic velocity. The observed offset however can be accommo- 
dated if there were a fairly modest misalignment between the 
disk axis and the jet. VLBI observations constrain the jet axis to 
be within 10° of the sky plane. For such an orientation misalign- 
ments between projected jet axis and disk axis will be similar 
to intrinsic misalignments. Given the arguments above that the 
H i absorption occurs at ~ 80 pc radius and give n the estimated 
central black hole mass of (Mbh = 2.5 • 10 9 Mm. Ta dhunter et all 
120031) the orbital velocity is v lot (r = 80pc) = 367 km s -1 . Given 
this orbital velocity a misalignment of only 21° is sufficient to 
explain the difference between the H i centroid and systemic ve- 
locities. 



4.2. Physical properties of the broad absorption component 
gas 

According to Mal onev et all (1 19961) the physical state of gas 
around an AGN is controlled by an ionisation parameter deter- 
mined from the ratio of the hard X-ray photon flux to local gas 
density. At a given radius r this ionisation parameter equals 

£ eff = L x ■ n x ■ r 2 ■ N 2 ™ (2) 

where Lx is the luminosity in > 2 keV X-rays and n is the gas 
number density. The final term takes account of the effects of 
X-ray photoelectric absorption where N22 is the the column den- 
sity in units of 10 22 crrT 2 (In Cygnus A X-ray observations give 
an estimate of N22 = 20 along the line of sig ht to the core. 
Youn g et al ., 2002). According to the model of Mal onev et al.l 
(1996, see their Fig. 3) for Cygnus A at r = 80 pc the gas 
fraction is predominantly atomic (> 90%) for a density range 
10 3 < n < 1.6 ■ 10 5 crrr 3 . The remaining < 10% of the gas is 
ionised at low densities (n = 10 3 cirT 3 ) and becomes molec- 
ular at higher densities (n w 10 5 crrr 3 ). Above densities of 
3.3 • 10 5 cm" 3 the gas is most ly molecular. 

The lMalonev et all 1 1996) model can be used to predict the 
disk Hi and free-free absorption opacity as a function of gas 
density; comparison with observations can then constrain the 
density. In order to predict opacities we need to convert gas vol- 
ume densities to column densities which requires estimates of 
the path length through the absorbing gas. Based on our estimate 
of the disk thickness of 20 pc and the maximum disk inclination 
angle with respect to the line of sight (Sect. 4.1) we estimate the 
geometric path length through the disk to be at least 31 pc. This 
is a strict lower limit and we adopt L eeom = 40 pc as a more 
likely value. Note that if the gas exists in clouds the effective 
path length through absorbing gas can be less than L geom but it 
cannot be larger. Using the above assumptions we find that up- 
per limits on free-free absorption toward the counter-jet do not 
give any useful gas density constraints. This is because, although 
the free-electron fraction increases at lower densities this is com- 
pensated for by a lower total gas column such that the free-free 
opacity continuum absorption at 21cm wavelength against the 
counter-jet is less than 0. 1 over all reasonable densities. 

For modelling Hi ab sorption the relevant parameters given 
bv lMalonev et all (1 19961) are the predicted gas temperature and 
atomic fraction versus effective ionisation parameter £ e g . At the 
radius of the absorbing Hi assuming densities high enough to 
give a predominately atomic (and not ionised) column, the in- 
crease of r sp i n due to radiative excitation fr om the radio core is 
negligible (see e.g. Ba hcall & EkersL 1 19691) . Given this we as- 
sume the atomic gas is thermalised and that its spin temperature 
(r sp in) equals the gas temperature (T). The predicted Hi line 
width in this cases is proportional to NmT~ l where Nm is the 
Hi column density and T the temperature. Both of the above 
quantities can be calculated as functions of the local density n. 
The former quantity equals the product of the total column den- 
sity N - n- L geom (for a uniform filled column) multiplied by the 
Hi abundance. This abundance is a function of f e jf (calculated 
using Fig. 3 in IMalonev et all 1 19961) which in turn is a function 
of n via eq. 2 (where N22 - N/10 22 ). Lik ewise T can be derive d 
as a function of «, again using Fig. 3 in (Malon ev et all ll 996). 
Combining the Nhi and T dependencies on n together we can 
predict A^iL -1 and hence the expected Hi absorption linewidth 
versus local density and then compare to observations. For den- 
sities n > 10 4 crrr 3 the predicted Hi absorption is larger than 
observed, this can however easily be reconciled with observa- 
tions if the Hi absorption is concentrated in clouds (see below). 
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At densities n < 10 4 cm 4 however there are no solutions that 
can provide enough Hi absorption to match what we observe, 
hence we can set a minimum density of 10 4 cirT 3 for our H i ab- 
sorbing gas. 

For the gas density range 1.0 ■ 10 4 crrT 3 < n < 3.2 • 10 4 cm -3 
there exist pure atomic phase gas solutions which fit the H i ob- 
servations. These solutions require the absorption to occur in 
clouds such that only a fraction /Los of a typical LOS passes 
through clouds. At the higher densities within this range the 
requirements on total gas column density are reduced because 
gas/spin temperatures decline making hydrogen atoms more ef- 
ficient at absorption. At the boundary in density between hav- 
ing predominately atomic or molecular phase clouds (at density 
4.6 ■ 10 4 cm 4 ) the spin temperature is ~ 114 K the cloud line 
of sight filling factor is / L0S = 0.032 and N - 1.82 ■ 10 23 crrT 2 
which is comparable to the X-ray absorption column along the 
line of sight to the central engine. At higher densities the clouds 
become predominately molecular but the observed H i absorp- 
tion can still be fitted by similar total gas column densities. The 
reason is that although the atomic fraction declines rapidly with 
increasing density the temperature also declines, this increases 
the efficiency of absorption per hydrogen atom which almost ex- 
actly compensates for the reduced abundance. 

In summary we cannot on the basis of our Hi observa- 
tions alone distinguish between models for the Hi absorbing 
gas phase which are primarily atomic or molecular. We can 
however set at minimum density (n > 10 4 cm 4 ) and a range 
for the total gas column density through the H i absorbing disk 
(10 23 cm 4 < N < 10 24 cm 4 ). Estimating total gas masses is 
complicated for those solutions invoking clouds in that these so- 
lutions only fix the line of sight filling factor, converting to vol- 
ume filling factors requires knowing the cloud size, on which we 
have no constraint. For the minimum density solution where the 
gas column is continuous we can make a rough estimate of disk 
gas mass of M easc iisk = 1.0 ■ 10 8 M Q within the inner 80 pc ra- 
dius of the disk. For higher density solutions, both total required 
column densities and cloud filling factors decrease so that to- 
tal gas mass requirements are significantly less. In all cases the 
total disk mass is much less than that of the c entral black hole 
(M BH = 2.5 ■ 10 9 Me, iTadhunter et al.L [2003b which therefore 
dominates the kinematics in the inner part of the nucleus that we 
observe. 



4.3. Constraints on circumnuclear torus properties and 
relation to the H i disk 

It is interesting that the estimated total gas column density 
through the H i absorbing gas (see Section 4.2) is comparable 
to that estimated by X-ra y photoelectric ab sorption toward the 
core of 2 ■ 10 23 cm 4 (see lYoung et al.L [20021) . Is it possible that 
most of this column density and most of the material that hides 
the quasar nucleus occurs as far out as r — 80 pc? An obvious 
problem with such a scenario is the absence of broad H i absorp- 
tion along the direct line of sight to the radio core (peak opacity 
is To = 0.016 ± 0,006, compared to > 1 on the counter-jet, see 
Sect. 3.3 and Fig. [2). A change in gas physical state with scale- 
height in the disk while keeping almost constant column density 
seems unlikely. If the high scale-height gas covering the core 
were still predominately atomic or molecular we find no den- 
sity solutions that fit the low limit on H i absorption seen. On the 
other hand if the high scale height gas were wholly ionised at this 
column density it would strongly free-free absorb the radio core. 
Although there are signs of free-free absorption at r < 20 pc 



(iKrichbaum et all 0.998) along the counter-jet this absorption 
does not cover the core position itself. 

It seems likely instead that most of the X-rays absorbing col- 
umn in Cygnus A and other hidden quasars occurs on scales 
<sc 100 pc in a compact circumnuclear torus. The inner radius 
of such a torus is set by th e dust sublima t ion rad ius, = 0.4 • 
(Lboi • lO -45 erg" 1 s) 5 pc dNenkova et all l2008bl) where L bol is 
the AGN bolometric lumi nosity. For Cygnus A Lb oi is estimated 
to be 1.5 • 10 45 erg s 4 ( W hvsong & A ntonucci. 2004) giving 
rg — 0.5 pc. It should be noted that the IWhvsong & Antonuccil 
(2004) mid-IR estimate of the bolometric luminosity is signif- 
icantly lower (by factors 5 - 20) than estimates based on th e 
unabsorbed hard X-ray luminosity (see ITadhunter et ail 12003). 
possibly because of uncorrected torus extinction at mid-IR wave- 
lengths (C. Tadhunter, private communication). However, an in- 
crease by a factor of 20 in bolometric luminosity results in the 
sublimation radius increasing only to r& = 2.2 pc, still signifi- 
cantly smaller than the radius of the H i structure. The outer ra- 
dius (r out ) is probably not as cl eanly defined as in some early de- 
pictions of doughnut-like tori dPadovani & UrrvLll992h and may 
even be continuous with larger circumnuclear disks (see refer- 
ences in Sect. 1). Consistent with these more general geometries 
the term 'torus' can be thought of, if one prefers, as an acronym 
for 'Thick/ toro i dal Obscurati on Required by Unified Schemes' 
dConwavl 1 1 9991; lEhtzuH 120081) rather than referring to a specific 
fixed geomet ry. 

Recently iPrivonl (l2009i) has fitted the observed Spectral 
Energy Distribution (SED) of the Cygnus A nucleus with a com- 
bined synchrotron jet, starburst an d torus model. The toru s fit- 
ting assumed the clumpy model of iNenkova et all d2008allbh . A 
wide range of solutions were obtained depending on assump- 
tions about the disk inclination and synchrotron jet spectrum. 
Most plausible solutions however had cloud number density per 
unit volume declining as r~ q with q — 1 and ratios of outer 
to inner radius, Y = 30 (giving r out = 18 pc for Cygnus A). 
Satisfyingly, the predicted total column density matches the X- 
ray absorption estimate within the errors. Despite the relatively 
large formal outer radius the concentration of clumps toward the 
inner edge means that most of the column density is concentrated 
at small radii. For a radial exponent of q — 1 half of the total col- 
umn density occurs within yfY ■ r m corresponding to 2.8 pc, i.e. 
on a much smaller scale than our observed H i absorption. 

The above results beg the question of relationship between 
the H i absorbing disk and the circumnuclear torus in Cy gnus A . 
The m ass in the torus is very low (6 ■ 10 5 M , following lElitzurl 
120081) and the lifetime of large scale-height clouds within it 
is sho rt (due to intercloud collisions, see e.g. iKrolik & Lepp. 
119891) . Hence, both to fuel the quasar-like nucleus at the ex- 
pected rate of ~ 1 M©yr and to replenish the torus clouds a 
much larger reservoir of material is required, which could reside 
in the H i absorbing disk. The mass in such a disk (see Section 
4.2) is sufficient to power the source for 10 7 to 10 s yr and may 
form part of a 'feeding structure' which funnels gas from kilo- 
parsec scales into the central black hole. A feeding connection 
between disk and torus seems likely but the mechanisms of gas 
transport between the two scales is unclear. 

Pertinent to the feeding question is whether the torus and the 
Hi absorbing disk are a continuous structure with one gradu- 
ally melding into the other, or if they are dynamically distinct 
and contain gas with very different physical conditions. The 
former possibility is motivated by the fitted torus outer radius 
'"out = 18 pc which is of the same order of magnitude as the 
radius of the Hi absorption (r - 80 pc). Furthermore the torus 
outer radius determined from SED fitting is not well constrained 
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by the data and may extend beyond this limit. However, in such 
a continuous model for clouds with approximately fixed density 
the clouds would likely be in an increasingly molecular state 
over a projected radius from 1 to 70 pc as the excitation pa- 
rameter decreased. It is hard to see how the rapid gradient in 
H i absorption seen along the counter-jet could be reproduced in 
such a case. Furthermore the observa tion of free-free absorption 
at intermediate positions (r < 20 pc, iRrichbaum et all 119981) is 
inconsistent with such a continuous model. It seems more likely 
that the torus — although ultimately fed from gas in the disk 
— is a distinct structure, perhaps generated in a n accr etion disk 
wind as in the model of lElitzur & Shlosmanl (2006). Clearly, 
more work is required to understand how the observed com- 
ponents in Cygnus A and other radio galaxies of 100 pc scale 
H i absorbing disks, inner ionised gas and circumnuclear tori are 
connected in self-consistent structures which both obscure and 
feed the central engine. 

An additional remark concerning the disk/torus inter- 
relationship is that although in Cygnus A the H i absorbing disk 
likely does not contribute to the total absorbing column toward 
the central engine (as estimated from the X-ray observations) it 
could have done so if we had observed Cygnus A from a differ- 
ent direction. 

Further progress on understanding the circumnuclear gas 
environment in Cygnus A probably requires the reliable de- 
tection of molecular absorption and its imaging with VLBI. 
Single dish searches at commonly observed molecular transi- 
tions such as the ground and higher rotational trans i tions of CO 
dBarvainis & Antonuccil Il994t ISalome & Combesl 120031) have 
so far yielded only upper limits. In contrast a tentative detec- 
tion of CO+ in abso rption using the IRAM 30m was reported by 
iFuente et alj (12000) with a centroid velocity and FWHM velocity 
width (170 km s ) very similar to that of our broad Hi absorp- 
tion. This detection has yet to be confirmed interferometrically 
(A. Fuente, private communication). Given the results in this pa- 
per, despite the similarity in spectral shape, it seems unlikely that 
these observations are probing the same gas column as seen in 
H i absorption. The counter-jet which provides the background 
continuum against which the atomic hydrogen absorption is 
seen will have negligible flux density at millimetre wavelengths. 
Instead the line of sight probed at 3mm wavelength will probably 
lie at projected distances < 1 pc from the central engine. While 
this is comparable to the scales on which the circumnuclear torus 
has its largest column density (~ 3 pc) we would expect that the 
velocity dispersion of clouds in such a geometrically thick torus 
would be comparable to the orbital velocity at this radius (i.e. 
1900 km s _1 ), which is much broader than is observed. Similar 
considerations apply to the ten tative VLBA observation s of ex- 
ited OH absorption reported bv llmpellizzeri et all d2006l) at pro- 
jected radii < 3 pc which only have observed width ~ 100 km s _1 . 

4.4. The narrow absorption component 

While we argue that the broad absorption component is caused 
by gas rotating around the central black hole (Sect. 4. 1) the nar- 
row (< 30 km s see Sect. 3) absorption gas is likely to have 
a different origin. This narrow velocity component is signifi- 
cantly redshifted (~ 186 km s _1 ) with respect to the optical 
systemic velocity implying foreground gas moving inward to- 
ward the nucleus. This component is detected over the whole 
continuum region (i.e. jet, core and counterjet) where the obser- 
vations are sensitive enough to detect an opacity of ~ 0.1 (see 
Sect. 3.2 and Fig. [2]). A very similar opacity and centroid ve- 
locity is seen at the core position and one effective beam away 



on the jet side. The apparent increase in opacity at the velocity 
of the narrow absorption that is seen on the counter-jet side in 
Fig. [2] can be explained entirely by contamination by the broad 
velocity wings of the broad velocity absorption. While we do 
not have many independent samples of the narrow velocity com- 
ponent opacity versus position those we do have are consistent 
with a constant narrow component opacity of 0. 1 over the whole 
VLBI radio source. Although the foreground distance from the 
narrow H i absorption to the VLBI continuum source is not well 
constrained we suspect that it must be > 100 pc given its narrow 
velocity width and the position stability of its velocity centroid; 
if the gas was closer one would expect tidal forces to widen the 
velocity width and increase velocity centroid variations. 

The structure and origin of the narrow H i absorption is un- 
clear. Physically it could be related to a minor merging even t, 
detected 400 pc South- West of the core dCanalizo et aU l2003h . 
As a result of this interaction the narrow H i component could 
arise in a tidal tail of gas that is movin g towards the nucleus . 
According to galaxy merger simulations dBournaud et all [2005) 
part of the progenitor's gas that gets expelled in tidal tails during 
a merging event will eventually fall back on the disk. The narrow 
component could also be connected to the giant infalling molec- 
ular cloud located 1 .35 kpc to the North- West of the nucleus 
dBellamv & Tadhunteill20 04). However, the projected distances 
from the core are rather different such that a direct connection 
is not obvious. Both of the above gas components are relatively 
large (> 100 pc) and thus would be consistent with having fairly 
constant Hi opacity over the whole VLBI structure. 

Further progress with constraining the size and origin of the 
narrow H i absorption system might be made by making future 
sensitive (and high spectral resolution) e-MERLIN or EVLA ob- 
servations to try to trace the narrow H i absorption further along 
the jet. 

5. Summary 

We have presented VLBA Hi absorption data of the core re- 
gion of Cygnus A. Hi absorption is detected over a linear 
scale of 95 pc, but is seen mainly along the counter-jet and on 
the nucleus. The integrated spectrum can be well-fitted by two 
Gaussian profiles suggesting abroad (FWHM= 231+21 km s -1 ) 
and a narrow (< 30 km s ) component. Modelling the data 
shows that the broad absorption occurs only against the counter- 
jet and not against the jet. Against the unresolved core the opac- 
ity is very low. The broad velocity component can be explained 
by a circumnuclear disk which has its highest opacity 45 mas 
away (in projection) from the black hole. The narrow velocity 
component could be explained as an infalling tidal tail, presum- 
ably left from a past minor merging event. 

The radius at which the broad H i absorption occurs can be 
constrained with the help of limits on the orientation of the disk. 
We find a relatively narrow range of possible parameters, result- 
ing in an estimated radius of ~ 80 pc for the peak opacity, a disk 
scale height of about 20 pc and hence an opening angle of 14°. 
The offset between the centroid of the broad velocity compo- 
nent and the mean systemic velocity can be explained by a tilt of 
21° of the disk axis compared to the jet axis. 

Based on the derived geometry of the circumnuclear disk we 
derived physical properties of the H i absorbing gas. We find the 
minimum gas density to be n > 10 4 cm~ 3 with a spin temperature 
^spin < 740 K. With the H i observations alone we cannot distin- 
guish between models for the H i absorbing gas phase which are 
primarily atomic or molecular. We can however set a range for 
the total gas column density through the H i absorbing disk of 
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10 23 cm -2 < N < 10 24 cm" 2 . An upper limit on the gas mass 
within a radius of 80 pc is M gas _disk = 10 8 M , which is a factor 
25 les s than the black hole mass estimated by Tad hunter et alj 
(120031) . 

The circumnuclear torus in Cygnus A has an estimated fidu- 
cial radius ~ 3 pc, which is a much smaller scale than our ob- 
served Hi absorption. The estimated mass in the torus clouds 
(6T0 5 M ) is too low to power (alone) the source for 10 7 — 10 8 yr. 
In contrast the H i disk has enough mass to feed the AGN and re- 
plenish the torus clouds. 

Higher sensitivity, broader band observations are needed 
to study the properties of the atomic circumnuclear disk in 
more detail, to search for a connection with the torus and 
to investigate whether the torus contains a very broad (2000 
km s _1 wide) H i absorption component. Higher spectral resolu- 
tion e-MERLIN or EVLA observations of the narrow absorption 
component are needed to try to trace absorption further out along 
the jet to constrain its physical size. 
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